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Detailed simulations of the gamma-ray sky played an important role before the launch of Fermi. Pre-launch 
simulation campaigns aided the development of new analysis tools and the assessment of the expected Fermi- 
LAT science performance. After four years of the Fermi mission in orbit, improved sky simulations can still 
be very useful to support the activities of Fermi, to refine our understanding of the scientific performance and 
to explore potential new science cases. We have used the LAT observations collected so far to develop new, 
highly-detailed simulations of the sky as seen by LAT. These simulations include the isotropic and Galactic 
diffuse emission, as well as the list of point and extended sources detected by the LAT. Furthermore, our 
sky simulations account for transient and variable sources, including Gamma Ray Bursts, pulsars and flaring 
blazars. We report on the latest developments of this new sky simulation work, describing its potential benefits 
for improving the science of Fermi today and in the near future. 



1. INTRODUCTION 

Before the launch of Fermi in June 2008, the 
Fermi Large Area Telescope (LAT) Collaboration 
conducted extensive sky simulation campaigns. Those 
simulations were very important for the development 
of the analysis tools and the assessment of the science 
performance of the LAT (jAtwood et all [2009j j V After 
4 years of Fermi operations, sky simulations are still 
very valuable to study the science potential of the 
LAT, because they allow the possibility of performing 
specific analysis under controlled conditions. The 
main areas where sky simulations are useful today 
are: 

-Catalog-related projects, including studies on sensi- 
tivity thresholds or false detection rate; 
-Instrument-related studies and validations; 
-Developments and improvements of analysis tech- 
niques, such as pulsar blind search algorithms; 
-Developments of new analysis tools; 
-Population-related studies (e.g. contribution of faint 
sources to the diffuse emission); 

-Evaluation of alternative observing strategies, in- 
cluding dedicated pointed observations; 



2. THE LAT OBSERVATION SIMULATOR 

The fundamental tool for creating simulations of the 
7-ray sky as observed by Fermi-LAT is the Observa- 
tion Simulator (gtobssim), which is part of the Fermi 
Science Tools 1 . Gtobssim creates a set of output simu- 
lated data products (including photon events and ob- 
servation intervals), which are formatted into FITS 
files, exactly as the high-level real LAT data available 



on the Fermi Science Support Center web page 2 . In 
order to produce these simulations of detected pho- 
tons, gtobssim needs three classes of inputs. First of 
all, it is fundamental to specify a sky model includ- 
ing the 7-ray sources to be simulated. The details 
about each simulated source (e.g. sky position, flux, 
spectrum) are described by XML files, although some 
sources, like pulsars or Gamma Ray Bursts require ad- 
ditional ASCII files, containing for example the light 
curve templates. Furthermore, gtobssim needs the 
LAT Instrument Response Functions (IRFs), which 
include the Point Spread Function, the Effective Area 
and the Energy Dispersion. A parameterization of the 
IRFs as a function of the energy and incidence angle of 
the incoming photon is thus provided to gtobssim in 
FITS format. The last piece of information needed by 
gtobssim is the pointing history of the Fermi satellite, 
which is crucial to determine the portion of the sky ob- 
served by the LAT at any time. Gtobssim can accept 
as input the real pointing history of Fermi in order 
to create a simulated replica of the real observation 
history. If no pointing history is provided, gtobssim 
uses a standard rocking profile and produces a FITS 
file, which is formatted as the real pointing history of 
the LAT. 

The distribution of the simulated photons as seen 
by the LAT is obtained by folding the source models 
with the IRFs and taking into account the pointing 
history. The output products, both created in FITS 
format, are an event file (FT1) containing the list of 
simulated detected photons with their characteristics 
(including the reconstructed energy, position, time of 
arrival). In case no pointing history has been provided 
by the user, a simulated pointing history file (FT2), 
based on the standard rocking profile, is also gener- 
ated. 



2 http: / / fcrmi . gsfc. nasa.gov / ssc / 
x http: / /fermi. gsfc. nasa.gov/ssc/data/analysis/scitools/overview. html 
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Figure 1: gtobssim-based simulation of a 4-year LAT observation of the full sky. The map shows only 7 rays above 1 
GeV that converted in the front part of the LAT tracker. The simulation is based on the real LAT pointing history of 
the LAT and the P7V6_SOURCE Instrument Response Functions. 



It is important to note that gtobssim does not take 
into account the detailed propagation of the particles 
in the LAT (this is done using a tool called Gleam 
based on th e GEANT4 particle interaction codes, as 
described in lAtwood et all [2009} . 



3. THE SKY SIMULATION 
INFRASTRUCTURE 

Simulating a large LAT observation of the full sky 
requires large computational resources, in particular 
when the number of sources is large and the observa- 
tion timeline covers few years. Therefore, in order to 
support massive sky simulations, we developed a ded- 
icated t ask running on the SLAC Pipeline I I fram e- 
work (jAtwood et alj [2009) . IZimmer et ail [2012j |L 
This task runs parallel jobs, each one simulating a 
small observation window (e.g. 1 week) and running 
post-processing operations, including the creation of 
Good Time Interval (GTIs) and the generation of live- 
timecube files. 

We also developed a set of utilities to manipulate the 
input sky models or to post-process the simulations, 
such as merging the output files, validating the simu- 
lation and creating the relevant documentation. 



4. SIMULATING THE 4-YEAR SKY OF 
FERMI-LAT 



dataset to compare with real data. We included in 
the sky model many different classes of sources rep- 
resenting the 7-ray emitters observed by the LAT. In 
particular, this sky model included the Galactic and 
the Isotropic diffuse emissio ns 3 , the emission from th e 
quiet Sun 4 and the Moon (jAbdo et al.l |201ll |2012| ) 
and the s ources in the Seco nd Fermi Source Cata- 
log f2FGL lNolan et al.1 [2012j ) i ncluding the extende d 
sources observed by the LAT (jLande et alj [2012) ). 
Furthermore, gtobssim is very flexible and allows the 
possibility of simulating also variable source s, includ- 
ing Gamma Ray B ursts ([Abdo et al.l 120091 1. blazars 
Nolan et al.l [2012j |) and 7-ray pulsars ()Razzano et al.l 
20091 ]). The contribution of unresolved Galactic and 
extragalactic sources produced by population s ynthe - 
sis cod e s 5 c an also be added ( Abdo et al.l 2010], 
IStrond [2007( h Figure 2 shows a comparison be- 
tween simulated and real data using a dataset of 6 
months. The simulation has been created using the 
P7V6_SOURCE IRFs and the real pointing history of 
the LAT. The comparison of the distributions of pho- 
ton characteristics (e.g. the energy) shows that the 
simulations agree with the data within few percents. 
This small discrepancy can be due to the absence of 
the Earth limb model in the simulation and does not 
affect the overall quality of the sky simulations. We 
therefore plan to use this sky simulation infrastructure 
to address more specific studies in the perspective of 
better understanding the science capabilities of the 
LAT in the next years of operations. 



We used the infrastructure described in Section 3 to 
create a realistic simulation of the 7-ray sky as seen 
by the LAT during its first four years of operations. 
This simulation (Fig 1) has the purpose of validat- 
ing the entire infrastructure and providing a reliable 



3 http: / /fermi. gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels. html 
4 http: / / sourceforge.net / projects / stellarics 

5 http:/ /www. mpe.mpg.de/ aws/talks/fermi_symp_strong_poster_ 
sourcepops.pdf 
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Figure 2: Comparison between the energy distributions of real LAT data and gtobssim simulation. This comparison is 
based on a 6-month simulation. The datasets agree within few percents and the discrepancy is likely due to the absence 
of the Earth limb emission in the simulations. 
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